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On the Most General Form of the Compatibility 
Equations and the Conditions of Integrability of 
Strain Rate and Strain 

Edmund H. Brown 

The most general forms of the equations of compatibility are derived by a very simple 
analysis in terms of the material strain rate tensor and the rate of change of the Riemann- 
Christoffel curvature tensor during the deformations. It is shown that recent statements 
on the significance of the compatibility equations, essentially that their satisfaction is 
equivalent to the condition that the space be locally Euclidean, are misinterpretations based 
on restricted forms of the equations. In addition, a new form of the conditions of inte- 
grability of strain is derived in terms of the rate of change of the curvature tensor. 

1. Definitions and Notations 

Because of the great variety of usage of authors in the field, a few of the definitions and 
notations employed here will be discussed briefly. However, a general knowledge of tensor 
analysis and modern elasticity is assumed. For the most general, and at the same time most 
readable introduction to tensor analysis, including applications to elasticity, the reader is 
referred to Brillouin [l]; 1 for an excellent and thorough review of modern elasticity through 
1953, despite a frequently cumbersome notation, the paper of Truesdell [2] is highly recom- 
mended; for more recent general discussions, the treatise of Green and Zerna [3], or the article 
of Doyle and Ericksen [4] should be consulted. 

To introduce the strain tensor, we imagine a complete //-dimensional material continuum 
G which, with an arbitrarily chosen coordinate system x\ is defined by the form of the metric 
tensor gij=gij(x k ). To show the covariant character of some of the quantities involved it 
will occasionally be necessary to transform from one coordinate system to another. Following 
Craig [5] we shall indicate quantities in one coordinate system by use of indices from the 
group i, j, k, h, I, and quantities in a second coordinate system by indices from the group 
Pi 2> r j s i t- Thus, such a transformation will be defined by an analytic relation of the form 
x r =x r (x k ). Further, for convenience, the components of the transformation matrices dx k l()x r 
or bx r /bx k will be denoted by Xr = dx k /dx r or X r k = c)x r ldx k . Similarly, &x k /i)x 8 dx r will be denoted 
by X k TS . During such a coordinate transformation, then, in one physical state, the components 
of the metric tensor will transform according to the tensor law, (j rs = g ij X i T X i s , the Christoffel 
symbols, or connections, will transform according to the Christoffel law, T* jk =T* st X s j X t k X i T -{- 
X s jk X\, and the components of the Riemann-Christoffel curvature tensor by the tensor law, 
R P q\rs = Rij\khX i p X j q X k r X h s . Note that we use a short vertical stroke between indices to emphasize 
various symmetry properties. In addition, following Green and Zerna, we use a full vertical 
stroke behind a quantity, as in the identity g ti \ k ^s0 9 to indicate the covariant derivative. 

A material coordinate system is defined b\ r assigning a value x % to each particle of the 
material continuum, which it retains during any subsequent deformation. We shall assume 
that the deformation, itself, depends on a single scalar parameter /x, independent of the 
coordinates x k , which in dynamic problems can be identified with the time t. Since the mate- 
rial coordinates of a particle do not change during a deformation, we must have £*(m 2 )=#*(aii). 
However, during the progress of the deformation the intrinsic quantities associated with 
the space, such as the metric tensor g ijy must change, assuming a continuous set of values 
9ij(pi)> Quite), .... When we wish to indicate several stages during a deformation, the 
appropriate values of p u /x 2 , . . . will be used, but more often only two stages, an "initial" 
state and a "final" state, are of interest, so that, for instance, we can write for the "initial" 
metric tensor, g iJ} and for the "final" metric tensor g iS . In such cases, we shall make a general 
practice of indicating a final state by placing bars over quantities which in the initial state 
are unbarred. 

1 Figures in brackets indicate the literature references at the end of this paper. 
443157—57 5 421 



> 



At this time, it is well to note that, parallel to the "material" representation of quantities, 
a "spatial" treatment can be used: A complete ^-dimensional geometric, or spatial, continuum 
H can be specified, using an arbitrarily chosen fixed spatial coordinate system y\ by the form 
of the spatial metric tensor, h ij =h ij (y k ). We imagine that the deformation of G takes place 
in H, that is, that there is a one-to-one correspondence between the "particles" of G and the 
"points" of H which, however, varies in form with ax. Thus, there is a "mapping" of G on H 
by y i =y i (x k ,fx) ) whose variation with /x represents the deformation of the material continuum. 
Since during the deformation, that is, during the variation of /x, the particles of (7 move through 
the points of H, 2/ t "(x fc ,jU 2 ) 7^(2*7x1), but at each given point y\ spatial quantities such as h tj and 
the spatial Riemann-Christoffel curvature tensor S^\ kh remain constant. There is nothing 
of course to prevent one assuming that the space H may change with parameters other than /x, 
a device which could possibly be of some use in discussions of plasticity; but for simplicity 
we shall avoid this here. In any case, a necessary condition that G remain in H is 

K W»-bp*\rs bx i dx i dx * dx n 

at each value of /x. 

Quantities such as strain can be related to either the material or spatial systems, whereas 
the displacement u* is essentially spatial. Since we intend to discuss strain only, for con- 
venience we will use chiefly a material representation and make little further reference to the 
spatial treatment. We should remark, however, that one of the principal applications of the 
results for non-Euclidean continua is in cases where (unlike our definitions above) the dimen- 
sion of G is m, the dimension of H is n, and m<w, since they are useful in reduced problems 
involving shells and curved rods. These problems, however, require the introduction of 
n—m other mth-order quantities ("changes in curvature", "changes in torsion", etc.) in 
addition to the strain, so that we shall not consider such possibilities here. 

2. Introduction of Strain and Derivation of Results 

We shall take as a measure of finite strain the quantity 

a2) e ij (x k ) = %[g ij (x k , l jL 2 )-g ij (x k ^ l )] ) (1) 

or, if we wish to think only of some particular initial and final states, 

etj(x*) = %[g t i(x k )-9ij(x k )], (2) 

since — even though it is less used in the modern theory of finite elasticity — it reduces to the 
familiar strain tensor of "classical" theory, such as is found in Love [6], in the infinitesimal case. 
According to the above definition, strain depends merely on the initial and final states, and is 
independent of the way in which the deformation took place. Such a definition is too general for 
strain to be uniquely related to stress in more complex cases involving permanent set or plastic 
behavior. Thus, we are led to the introduction of the rate of strain tensor e tj defined by 

. __&</_ 1 i>g tj , . 

We note that, since differentiation with respect to the scalar parameter /x does not affect the 
character of g ij} e w like g tJ is a tensor in the x k coordinate system for a given value of /x. Sim- 
ilarly, since both g tj and g Xi are tensors, the finite strain a2) e tj is a tensor. However, though 
x k =x k for all /x> the identities (based on the hypothesis of a Kiemann space, that is, that a 
unique measure exists) g i:j \i = and gfy|*=0jdo not allow us to conclude either that g tj \i^0 } or 
7/^.| fc ==0, since g tj is not the metric tensor of (?, nor "g tj the metric tensor of G. In G, g^ is merely 
a tensor field and not the metric tensor; similar statements can be made about the character of 
g i} in G. 
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Now, considering e fj to be a tensor field in G, in accord with the usual requirements for 
raising and lowering indices, we must have 

iti=0*g*i u =-l%£! (4) 

and 

Just as the metric tensor g tj changes during the deformation, the Christoffcl connection 
T* jk based on g i} must change. We shall define the rate of change Ej k of r*# by 

However, since T* jk transforms, when the coordinates are changed at one physical state, 
by the Christoffel law, and none of the quantities x k , X k , X k s are functions of ju, 

E: st =^ = c ^XiX k t X\=EU k XlX k X\ ) (7) 

which shows that — unlike T* jk — the quantities Ej k , and thus En Jk =g ih E* jk , are tensors, a fact 
first pointed out by Bergmann [7]. 

Just as T\ n can be expressed, for a Riemann space, in terms of the metric tensor, 

' }k 2 g [_dx k ~ t ~dx j dx h J 
so the E[ jk can be written in terms of the strain rate e tj ; for 

^ = &r! iA; 1 dg ih ri>g hJ dg hk d^l 1 d fdg h) dg hk Dg jk ~] 
' jk b M 2 d/i L*te* &*' dx h J r 2 J d f jildx k ~ r dx j dx h J w 

or, since #^ is an analytic function of x l and /i, and cc* and ,u are independent, 

E jk - 2e T hljk +g |^— +— -— J. (9) 

Similarly, Ei\ jk can be shown to be given by 

^i*=ff^«=-2^r?„+[^+^-||], (io) 

or 

Ei\j k =— ^ 2eihT.j k . (11) 

Using the expression for the covariant derivative of the second order tensor e ih the three 
partial derivatives of e# appearing above may be written 

^r-J == ink\j J r^iJcT m hj~\-€hiT l tkj 

"">I5 ,=a ""^i»l»""^»r./»— %r* ftft . (12) 

423 



On adding these expressions, interchanging indices where symmetry permits, we obtain 

§+^-|j=^U+4*L-^| ft +2e w r' M . (13) 

Introducing this in the expressions for E\ jk and Ei\ jk in terms of the strain rate, and making 
use of the fact that g ih \ k = results in the more elegant forms, 

E i \ jk =e ij \ k +e ki \ j —e jk \ i (14) 

and 

&j*=eUU+ik i \ j -~e jk \ i . (15) 

We must, finally, consider how the Riemann-Christoffel curvature tensor R\j\ k - n based on 
g tj varies during the deformation. Denning dT^iA^^dK^i^/djii, this rate of change of curva- 
ture is given by 



ar 1 * „, , ^ &r'» br l , jk , . dr 



or 



==^+^^ (16) 

<o] j\kh = E) jh\k lii % .jk\hy (I') 



and similarly, 



<^il*a — 9il&.j\kh — ^ ^ilR.j\kh — £Li\jh\k Ei\jk\h) (18) 



both of which are especially elegant forms for the rate of change of the curvature tensor. 

Just as in the case of E t \ jk , we can express S'a^h in terms of the strain rate, and, in doing so, 
we obtain the equations of compatibility in their most general form. By substitution and 
simple algebra these identities may be written 



?ij\kh~ 



: ■(««!»*— iij\kh)+(fith\jk— tain)— (tjh\ik— e jk \ ih ). (19) 



It is well to examine the significance of these equations. There is no reason for us to 
assume that the strain rate is a tensor function (which would require that e tj satisfy certain 
conditions of integrability), for in the most general case, the e ti would be merely the com- 
ponents of a tensor field. However, the curvature tensor Rij\kh is an intrinsic property of the 
material continuum G and must be related, as stated previously, to the curvature tensor S t j\ kh 
of the spatial continuum H, in order that G remain in H. Thus, #ij\ kh must be related to S iMh 
by an equation of the form 

6.m-*. nr« fylzypwbtf d?J (20) 

We have assumed that H does not vary, that is, that the value of S tj \ kh at each point y k is 
constant. In particular, if S ism =0 } R ijm =0, and thus, of course <£ ij]kh must also vanish. If, 
however, S tj \ kh is not zero, # mh will be constant, or more strongly, zero, only if other special 
conditions hold. 

The conditions of integrability of the rate of strain tensor must now be investigated. 
Using the relation for ^m, and the fact that Rij\ kh is skew-symmetric in the indices of both 
pairs, we can write 

&ij\kh^ &ji\kh— ~2[&ilR. j\kh~\- &jlR. i\kh\- (21) 
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However, the equations 

VilR.j\kh~hejlft.i\kh = eij\kh ^ij\hk (22) 

are merely the conditions of integrability; the vanishing of either side of the above equation 
implies that if e tj is subjected to parallel transport around a closed curve, it will return to its 
initial value at the starting point, and thus, may be considered a tensor function of the coordi- 
nates x l rather than merely a tensor field. 

Thus, the rate of strain is integrable if any one of the following equations holds identically: 

£ij\kh &ij\hk == 

e i iR! J \kh+e 3l R! i \kh=0 (23) 

(oij\kh~r Sji |fcft= 0, 

the last of the equations being equivalent to the requirement that Stj\kh be skew-symmetric 
in its first two indices (it is already skew-symmetric in the second pair as a consequence of its 
definition). It should be noticed that these integrability conditions are distinct from the 
compatibility equations, and in particular, it is not necessary for Sij\m to vanish for e tj to be 
integrable. 

If, however, the strain rate is integrable, and we substitute the first of the above conditions 
in the compatibility equations, they take the simpler form 

<oij\kh = \&ih\jk ^ik\jh) \&jh\ik — ^jk\ih) • (24) 

In a great many practical cases, either the space is Euclidean, or some special relation 
may hold, such that Sa\m=^ ) which reduces the compatibility equations to the form 

6ih\jk 0ik\jh = 6jh\ik e jk\ih- (25) 

It is of some interest to see the forms which the last equations take when n<3. If n=l, 
all indices must be identical and we have, merely, = 0. If n=2, only two indices can be 
distinct, resulting in only one independent equation 

e ii \kk+e k jc\ i i=2e ik \ ik) i^k, (26) 

indices not summed. 

If n=Z, three indices may be distinct, and two independent sets of, in all, six equations 
maybe derived: 

eii\kk+e hk \ii=2e ik \a i^k, (27) 

indices not summed (3 equations), and 

««|/»+«i*|«=««|«+«<*|*yj i^j^k, (28) 

indices not summed (3 equations). 

The last two sets of equations, for n=3, are equivalent to those given by A. E. H. Love 
in Cartesian coordinates. Thus, our completely general form of the compatibility equations 
reduces to the more well-known form for the appropriate conditions. 

We note that, since /x is a scalar parameter independent of x\ derivatives of the infini- 
tesimal strain tensor de ij =e ij diJL may be written dei j \j c =e ij dfx\ k =e ij \ k diji. Similarly, derivatives 
of the finite strain &ij=S&iflV' may be written e id \ k -= f eijdix\ k ^= f e id \ k d\x. Thus, all the rela- 
tions and conditions derived above hold equally well for both the infinitesimal and finite strain 
tensors. 

In conclusion, we can state the equations of integrability and the equations of compati- 
bility are two independent sets of conditions. Satisfaction of the'! former implies that the 
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space to which the deforming medium is assigned is Euclidean. Satisfaction of the latter 
implies the medium remains in that assigned space. The fact that both must hold for the 
solution of actual problems in elasticity should not lead to considering them as identical, for 
each has its own physical significance. 
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